function [h1,h2] = plot(spectrum,varargin)
% if nargin == 1
%     varargin = {'bo','linewidth',1};
% end
% h = stem(spectrum.mode_id,spectrum.mode_value,varargin{:});
mode_value = spectrum.mode_value;
mode_id = spectrum.mode_id;
if(nargin==2)
    mode_id = varargin{1};
end
subplot(2,1,1);
h1 = stem(mode_id,mode_value(1,mode_id+1));
set(gca,'FontSize',18);
% xlabel('mode ID');
ylabel('\phi_1 (sin)');

subplot(2,1,2);
h2 = stem(mode_id,mode_value(2,mode_id+1));
set(gca,'FontSize',18);
xlabel('mode ID');
ylabel('\phi_2 (cos)');